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•^ Abstract: Bolometers are phonon mediated detectors used in particle physics experiments to 

• search for rare processes, such as neutrinoless double beta decay and dark matter interactions. 

O They feature an excellent energy resolution, which is a few keV over an energy range extending 

^ from a few keV up to several MeV. Nevertheless the resolution can be limited by the noise induced 

(~| by vibrations of the mechanical appai^atus. In arrays of bolometers part of this noise is con^elated 

, ^ among different detectors and can be removed using a multichannel decorrelation algorithm. In 

this paper we present a decorrelation method and its application to data from the CUORICINO 

►^ experiment, an array of 62 TeOi bolometers. 
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1. Introduction 

Bolometers are detectors in which the energy from particle interactions is converted to thermal 
energy and measured via their rise in temperature. They provide excellent energy resolution, 
though their response is slow compared to conventional detectors. These features make them a 
suitable choice for experiments searching for rare processes, such as neutrinoless double beta de- 
cay (OvDBD) and dark matter (DM) interactions. 

The CUORE experiment will search for OvDBD of '^°Te QI], 0] using an array of 988 TeOi 
bolometers of 750 g each. Operated at a temperature of about 10 mK, these detectors feature an en- 
ergy resolution of a few keV over their energy range, extending from a few keV up to several MeV. 
At low energies, below 100 keV, the resolution is of the order of 1 keV FWHM, while at 2528 keV, 
the OvDBD energy, is about 5 keVFWHM; this, together with the low background and the high 
mass of the experiment, determines the sensitivity to the OvDBD. CUORE could also search for 
DM interactions, provided that the energy thi^eshold is of few keV [2D- An experiment made of 
62 bolometers, CUORICINO \M, was operated at Laboratori Nazionali del Gran Sasso (LNGS) in 
Italy between 2003 and 2008, and proved the feasibility of the TeOi bolometric technique. CUORE 
is currently under construction, and will start taking data in three yeai^s. 

Vibrations of the cryogenic apparatus (keeping the system at low temperatures) induce a visible 
noise, which limits the energy resolution at low energies [Q] and the energy threshold [Q]. Since all 
bolometers are held by the same structure, part of the vibrational noise is expected to be correlated. 

In this paper we present a method to estimate the correlated noise among different bolometers, 
and a method to remove it. The application to data from CUORICINO shows that the correlated 
noise is visible and that it can be efficiently removed. 

The CUORE detector structure will be different from the CUORICINO one, and thus a differ- 
ent vibrational noise is expected. If it will be unfortunately large, the algorithms developed in this 
work could be used to improve the performances of CUORE in both OvDBD and DM searches. 
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2. Experimental setup 

CUORICINO and CUORE bolometers are composed of two main parts, a Te02 crystal and a neu- 
tron transmutation doped Germanium (NTD-Ge) thermistor [0, [7|]. The crystal is cube-shaped 
(5x5x5 cm"*) and held by Teflon supports in copper frames. The frames are connected to the mix- 
ing chamber of a dilution refrigerator, which keeps the system at a temperature of about 10 mK. 
The thermistor is glued to the crystal and acts as thermometer. When energy is released in the crys- 
tal, its temperature increases and changes the resistance of the thermistor. To read out the signal, 
the thermistor is biased with constant current, which is provided by a voltage generator and a load 
resistor in series with the thermistor. The resistance of the thermistor varies in time with the tem- 
perature, and the voltage across it is the bolometer signal. The value of the load resistor is chosen 
to be much higher than the thermistor, so that the voltage across the thermistor is proportional to 
its resistance. The wires that connect the thermistor to the electronics introduce a non-negligible 
parasitic capacitance. 

The typical response of CUORICINO and CUORE bolometers to particles impinging on the 
crystal is of order 100 /iV/MeV. The signal frequency bandwidth is — 10 Hz, while the noise 
components extend to higher frequencies. The signal is amplified, filtered with a 6-pole active 
Bessel filter with a cut-off frequency of 12 Hz and then acquired with an 18-bit ADC with a sam- 
pling frequency of 125 Hz. The gain of the amplifiers ranges from 500 to 10000 V/V, and is 
tuned for each bolometer to fit the signals in the ADC range, which is ±10.5 V. A typical sig- 
nal recorded by the ADC, produced by a 1461 keV 7 particle fully absorbed in a CUORICINO 
bolometer, is shown in the left panel of Fig. [I]. 
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Figure 1. Response of a CUORICINO bolometer (channel 13) . Left panel: signal generated by a 1461 keV 
7 particle recorded by the ADC. Right panel: noise power spectrum expressed in mV^/Hz, and energy 
spectrum of the signal scaled to fit in the noise range. 



The 62 bolometers of CUORICINO were arranged in a tower of 13 floors: 1 1 floors were com- 
posed by 4 5x5x5 cm^ bolometers each, 2 floors were composed by 9 small bolometers (3x3x6 cm^) 
each. The small bolometers were recycled from previous arrays, and will not be used in CUORE, 
which will only use 5x5x5 cm^ crystals arranged in 19 towers of 13 floors each (Fig. ^. The front- 
end electronics, which provide the bias, the load resistors and the amplifier, are placed outside the 
cryostat, at room temperature [Q]. In CUORICINO, 1/3 of the electronics channels used load re- 
sistors and pre-amplifiers located inside the cryostat, at 1 10 K. This was done to reduce the noise 
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Figure 2. Arrays of Te02 bolometers. The CUORICINO tower (left), composed by 62 bolometers, rep- 
resented inside a schema of the cryostat; a single CUORE tower (middle) and the entire CUORE array 
(right). 

from load resistors, but since no improvement was achieved, the cold electronics will not be used 
in CUORE. Three channels ( 1 cold, 2 warm) were lost after the initial CUORICINO cool-down, 
reducing the number of active chaimels to 59. 

Vibrations of the detector structure generate two types of noise: thermal and microphonic. The 
thermal noise is due to vibrations of the crystals that induce temperature fluctuations of the crystals 
themselves. The microphonic noise is due to vibrations of the wires that connect the thermistor to 
the cryostat socket. A typical noise power spectrum is shown in the right panel of Fig. |I]; the peaks 
in the figure are attributed to microphonism, while the continuum is attributed to crystal vibrations. 

The data acquision system connected to the ADC boards implements two kind of software 
triggers. The first trigger fires when a signal is detected on a bolometer, the second fires randomly 
in time to acquire noise waveforms. To estimate the noise correlation, the acquisition was pro- 
grammed to acquire simultaneously all the bolometers of CUORICINO when the signal or the 
noise triggers fired on a bolometer. The length of the acquisition window was set to M = 5 12 sam- 
ples on all channels, corresponding to 4.096 s. The data used in this paper amount to 2 days taken 
from the last month of CUORICINO operation. 

3. Noise correlation 

The noise power spectrum is estimated from a large number of waveforms acquired with the ran- 
dom trigger, removing those that, by chance, contain signals. The power spectrum Nj{(o) of each 
bolometer channel / is computed as: 



Ni{(0) =<ni{(0)-n*{(0) >, 



(3.1) 
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where <> denotes the average taken over a large number of waveforms and n{o)) is the discrete 
Fourier transform (DFT) of a noise waveform «,(?). To avoid DFT artifacts, each waveform is 
weighted before the DFT using a Welch windowing function: 
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In the above expression the parameter G is a normalization coefficient chosen to satisfy the condi- 
tion m- 
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where a, is the standard deviation of the noise in the time domain. 

The covariance between the noises on bolometers channels / and j is estimated as: 



Nij{(0)=<ni{0))-n*{(0) > 



(3.4) 



where the definitions are the same of Eq. pTl] . The magnitude of Nij{(o) is the covariance as usually 
intended for real variables, while the phase accounts for any possible time delay of the frequency 
CO observed in the two bolometers. The correlation is also a complex quantity, and is defined as: 
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As an example, we show the behaviour of p^ i(g)) in Fig. ^; the microfonic peaks and the low 
frequencies are highly correlated, while the intermediate frequencies are almost uncorrelated. The 
pattern in the figure is visible for many {i,j) pairs and indicates that the source of correlated noise 
has analogous origin in all bolometers. 

The average correlation over all frequencies can be computed as: 
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and is found to be around 60% for the {ij) pair (13,1). The average correlation matrix, including 
the correlation of all the CUORICINO channels, is shown in Fig. ^ The matrix is not exactly 
triangular because Eq. ^]^is not symmetric under the exchange of channels / and j, since the same 
frequency may contribute with a different weight to Ni{(o) and Nj{(o). The figure shows that the 
correlation is high, even tough there is no evident link with the spatial disposition of the bolometers. 
Bolometers close to each other, in fact, aie not necessarily the most coiTclated. 

We also tried to order the channels by their cabling group in the cryostat or by cold/warm 
electronics channels, but also in this case we did not see any link. Other variables linked to the noise 
correlation could be, for example, the path of the wires along the tower, the tightness of the crystals 
and the length of thermistor wires. These variables are not accessible now since CUORICINO has 
been disassembled. Moreover, while the understanding of the noise sources is important, in this 
paper we focus on the analysis algorithms to reduce the observed noise. 
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Figure 3. Power spectra of CUORICINO channels 13 and 1 and their correlation p(ft)). 



4. Noise decorrelation 

In Te02 bolometers the noise is purely additive to the signal, meaning that the waveform f{t) 
observed in presence of a signal of amplitude A can be expressed as: 



m=A-s{t)+n{t) 



(4.1) 



where s{t) is the signal shape. Since the noise on different bolometers is partially correlated, it is 
possible to remove part of n(f ) from f{t) using the waveforms from other bolometers in which only 
noise is present. In the following we develop a method to remove the correlated noise. 

Assuming that each frequency component of «,(&)) is normally distributed, the multidimen- 
sional probability distribution of the noises on all bolometers n{co) = {n\ {(o),n2{(o),. . .} can be 
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Figure 4. Noise correlation in CUORICINO. Channels are ordered as bolometers are ordered in the array. 
White lines indicate dead channels. 



written as: 



P[n{(o)] oc exp [ - -n^ {(0)N ^{(o)n{(o) 



(4.2) 



where A'^ is the covariance matrix whose elements are defined in Eq. [3.4| . The probability dis- 
tribution of the noise on bolometer /, «,(w), can be obtained integrating out the noise of other 
bolometers: 
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The above equation represents a Gaussian with mean 
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and variance 1 jN^^ ( ft) ) 



The decorrelated value of «,(g)) can be obtained for each waveform as: 



f («) = mico) + X ^^^M^'') ' (4.4) 

and its power spectrum is expected to be: 

Nf(co) = — -. . (4.5) 

A generic waveform /,(?) on bolometer /, containing noise or noise+signal, can be decorrelated 
using the noise from all other bolometers as: 

/;'^(«)=/,(«) + £j^ny(«)=A5,-(«) + nf(«). (4.6) 

Summarizing, once the covariance matrix in Eq. ^^ is estimated from the data, the waveforms 
on each bolometer are decon^elated using Eq. ^3| and the effect on the noise power spectrum can 
be predicted with Eq. 



5. Application to data 

The acquired data are split in two sets, one set is used to estimate the correlation matrix and the 
expected noise power spectra, the other is used to verify that the application of the decorrelation 
algorithm to waveforms produces results consistent with the expectations. Figure § (left) shows 
the original power spectrum of channel 13 (Eq. PTTD and the one expected decorrelating from all 
the other CUORICINO bolometers (Eq. P3| ). The power of low frequencies is reduced and the 
microphonic peaks are completely removed. In principle one could be satisfied and proceed to 
apply Eq. ^^ to the single waveforms. However, using all the bolometers to decorrelate every 
triggered waveform is expensive from the computational point of view, since it requires DFTs for 
each bolometer used. Moreover all the waveforms used to decorrelate should not contain signals. 
This requirement is often fulfilled in CUORICINO, where the overall counting rate was around 
0. 1 Hz, but not in CUORE, where the rate is expected to be lai^ger. For this reason we tested to 
see if the decorrelation is effective using a smaller number of bolometers. For each bolometer 
we selected the most con^elated bolometers and we computed the expected decorrelated power 
spectrum. The results obtained are equivalent to those obtained using all bolometers (Fig. ^ left). 
The number of bolometers used to decorrelate was set to 1 1, a number sufficiently high to maximize 
the decorrelation and sufficiently low to ensure good performances. 

The noise power spectrum obtained on waveforms decorrelated with Eq. ^3| is shown in Fig. ^ 
(right). The decorrelation is found to be very effective on channel 13: the noise a is reduced from 
2.0 mV to 1.2 mV (corresponding to 3.6 and 2.2 keV, respectively), an effect which is visible also 
in the time domain (Fig. ^. We show in Fig. [7| the energy resolution on all bolometers before and 
after the decorrelation: in a few bolometers the resolution remains the same, while in others it is 
reduced up to 50%. 
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Figure 5. Left: Noise power spectra of bolometer 13: original (solid black), expected decorrelated using 61 
bolometers (dashed black), and expected decorrelated using the 1 1 most correlated bolometers (solid red). 
Right: original noise power spectrum and power spectrum of waveforms decorrelated using Eq. |4-.6t The 
original power spectra in the figures differs slightly because different sets of data has been used to estimate 
the covariance matrix (left) and the effect of the decorrelation (right). 
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Figure 6. Left: original (solid black) and decorrelated (dashed red) noise waveform from channel 13. Right: 
signal generated by a 43 keV 7 particle. The decorrelation algorithm reduces the noise leaving the signal 
unmodified. 



6. Combination with the optimum filter 

In the data analysis of Te02 bolometers the signal amplitude is estimated using the optimum fil- 
ter [HO . \nW . a filter that significantly improves the energy resolution. This filter can be used when 
the shape of the signal and the noise power spectrum of each channel are known. The filter acts on 
a single channel and its transfer function is 



H{co) = h 



N{(o) 



WIM 



(6.1) 



where im is a parameter to adjust the delay of the filter. In this work, as in Ref. [2D, we chose 
it equal to the maximum position of Sj, so that the maximum of the filtered signal is aligned to 
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Figure 7. Noise resolution of all CUORICINO bolometers before and after the decorrelation. Values are 
expressed in keV. 



the maximum of the non-filtered one. /j is a normalization constant that leaves unmodified the 



N{m) 



amplitude of the signal: h = 

To combine the optimum filter with the decorrelation, it is sufficient to replace in Eq. ^TT] the 
original power spectrum N{co) with N''{(0) defined in Eq. 4^, and then process the decorrelated 



waveforms in Eq. PT^ The comparison of the filtered noise power spectrum with and without 
the decorrelation for two bolometers with opposite behavior (channels 13 and 38) is shown in 
Fig. |S[ In channel 13, while the decorrelated noise is lower than the original one, there is no big 
improvement in the filtered-decorrelated noise, because the correlation is high only on frequencies 
with low signal to noise ratio. In channel 38, on the other hand, the decoiTclation lowers the noise 
after the optimum filter. 

To summarize we compare the expected resolution of the entire array using both decorrelation 
and optimum filter with respect to optimum filter only (Fig. ^. The improvement is not significant 
as when we decorrelated the original waveforms (Fig. [7|), because in most cases the noise is not 
correlated in the signal frequency bandwidth. 



7. Conclusions 

In this paper we developed a method to remove the correlated noise between different detectors. 
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Figure 8. Comparison of original, decorrelated, optimum filtered and decorrelated/optimum filtered noises 
on bolometers 13 (left) and 38 (right). 
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Figure 9. Noise resolution of all CUORICINO bolometers after the optimum filter before and after the 
decorrelation. 



The apphcation to CUORICINO showed that the correlated noise, generated by vibrations of the 
detector structure, can be efficiently removed. However, when the decorrelation is combined with 
the optimum filter, the resolution does not improve significantly since the correlated noise lies at 
frequencies higher than the signal frequency bandwidth. CUORE, a TeOi bolometric array 20 times 
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larger than CUORICINO, will have a new detector structure that could induce different vibrational 
noise. Moreover, while CUORICINO was refrigerated using liquid helium, the CUORE cryostat 
will use pulse tubes, acoustic cryocoolers oscillating at low frequencies (around 2 Hz). If the noise 
will be unfortunately large in the signal band, the method we introduced will be a valid tool to 
improve the CUORE performances. 
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